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Abstract 

Using realistic low-energy electron models, derived from the first-principles electronic structure 
calculations, we investigate behavior of interatomic exchange interactions in Cr 02 , which is re¬ 
garded to be one of the canonical half-metalic (HF) ferromagnetic. For these purposes we employ 
the dynamical mean-field theory (DMFT), based on the exact diagonalization of the effective An¬ 
derson impurity Hamiltonian, which was further supplemented with the theory of infinitesimal spin 
rotations for the exchange interactions. In order to elucidate the relative roles played by static 
and dynamic electron correlations, we compare the obtained results with several static techniques, 
including the unrestricted Hartree-Fock (HF) approximation, static DMFT (corresponding to the 
infinite frequency limit for the self-energy), and optimized effective potential (OEP) method for 
treating the correlation interactions in the random-phase approximation. Our results demonstrate 
that the origin of the HM ferromagnetism in Cr02 is highly nontrivial. As far as the interactions in 
neighboring coordination spheres are concerned, HF and DMFT methods produce very similar re¬ 
sults, due to the partial cancelation of ferromagnetic (FM) double exchange and antiferromagnetic 
(AFM) superexchange contributions, which represent two leading terms in the (AS) -1 expansion 
for the exchange interactions (AS being the intraatomic exchange splitting between majority- and 
minority-spin states). Both contributions are weaker in the HF approximation due to, respectively, 
additional orbital polarization of the t 2 g states and neglect of dynamic correlations. The role of 
higher-order terms in the (AS) -1 expansion is twofold. On the one hand, they give rise to addi¬ 
tional FM contributions to the neighboring exchange interactions, which tend to stabilize the FM 
state. On the other hand, they produce AFM long-range interactions, which make the FM state 
unstable in the DMFT calculations for the minimal model, consisting of the t 2 g bands. Thus, the 
robust ferromagnetism in the t 2 g model, which can be easily obtained using static approximations, 
is fortuitous and this picture is largely revised at the level of more rigorous DMFT approach. We 
argue that the main ingredients, which are missing in the t 2 g model, are the direct exchange inter¬ 
actions and the magnetic polarization of the oxygen 2 p band. We evaluate these contributions in 
the local-spin-density approximation and argue that they play a very important role in stability of 
the FM ground state in Cr02- 


* SOLOVYEV.Igor@nims.go.jp 


2 



I. INTRODUCTION 


Cr02 provides a rare example of metallic ferromagnetism in stoichiometric oxides. It is 
widely used in magnetic recording and still considered as one of the best particulate ever 


invented for these purposes [l|, jj. Besides magnetorecording, chromium dioxide has at¬ 
tracted a considerable interest due to its half-metallic (HM) electronic structure, which was 

n 

predicted by first-principles calculations [3J. The HM electronic structure is such that the 
majority-spin electrons are metallic, whereas the minority-spin electrons are semiconducting 
4]. In CrC> 2 , such behavior has been supported by point-contact Andreev reflection mea¬ 
surements 5]. Because of its implication in various spin-dependent transport phenomena 
6], the half-metallicity is the very important property of magnetic substances, which is in¬ 
tensively studied today [7]. These studies typically include both fundamental and practical 
aspects. 

Needless to say that ferromagnetism is one of the key properties of CrC> 2 , which predeter¬ 
mines its popularity and importance in all the applications. Although the Curie temperature 
is not exceptionally high from the view point of practical applications (about 390 K, mean¬ 
ing that the magnetic properties are significantly deteriorated at room temperature) 0 , it 
is still sufficiently high in order to classify CrC >2 as “robust ferromagnet”. 

Because of its popularity, Cr0 2 is the well studied material, both theoretically and exper¬ 
imentally. There is a fair number of theoretical works, focusing on the analysis of structural, 
transport, optical, and electronic properties of CrC >2 jdl. Many of them are based on 

the first-principles electronic structure calculations. These works clarify many important 
aspects of the material properties of Cr02- However, despite its immanent importance in 
the field, the problem of interatomic magnetic interactions and stability of the ferromagnetic 
(FM) ground state is CrC >2 remains in the shadow. Particularly, why is CrC >2 ferromagnetic? 
What are the main microscopic mechanisms yielding the FM ground state in Cr02? From 
our point of view, these important questions remain largely unanswered and in the present 
work we will try to fill in this gap. 

The ferromagnetism in CrC >2 is typically ascribed to the double exchange (DE) mecha¬ 


nism 


10 


15 


, ll6j. which was originally introduced for magnetoresistive manganites 17-20]. 
This mechanism is governed by the large intraatomic exchange splitting (AS) between the 
majority (d) and minority (d) spin states, which penalizes the electron hoppings between 
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atoms with opposite directions of spins. In such situation, the FM state will be the most 
stable one because any deviation from the collincar FM alignment of spins will increase the 
kinetic energy of electrons. The DE picture is well justified for large-spin (S) systems. In 
manganites, where S — 2, it can be very useful for semi-quantitative analysis, and, in many 
cas es, p rovides a valuable insight in understanding the electronic and magnetic properties 


20 


21]. However, even in this case, the additional effects can play a very important role and 


substantially modify the canonical DE picture 


22], For instance, the well-known antiferro¬ 
magnetic (AFM) superexchange interaction 23], which is also important in manganites, is 
y a next-order effect in the (AS) -1 expansion for interatomic exchange interactions 


'on 

nal 

21, 

22 


. In Cr02, where S — 1, the interatomic exchange splitting is not particularly large 
and the DE picture can be even more problematic: namely, besides FM DE interactions 
between the nearest neighbors, one can expect other magnetic interactions (not necessarily 


21 


. Another important factor, 


the FM ones), which can alter the magnetic ground state 
which is not treated by the DE model, is the oxygen states |24M27]. 

Another disputable point is the role of electron correlations beyond the local-spin-density 
approximation (LSDA) and whether Cr02 should be regarded as a strongly-correlated ma¬ 
terial or not. On the one hand, LSDA and generalized gradient approximation (GGA) 
already provide a reasonable description for the structural, transport, and optical properties 
of Cr02 with only moderate manifestation of many-body effects [8, 9|, 111]. On the other 
hand, it was also suggested that electron correlations are essential for understanding results 
of photoemission, x-ray absorption, optical, and resistivity measurements |lj|. We are not 
aware of any investigation of the effect of electron correlations on the behavior of inter¬ 


atomic exchange interactions in Cr02- Basically, there is only one theoretical work 


which addresses the problem of interatomic exchange interactions in Cr 02 on the basis of 
first-principles GGA and LSDA+f/ calculations. However, both are static techniques and 
do not treat dynamic correlations. Moreover, the reliability of the LSDA+D approach suf¬ 
fers from the use of adjustable parameters and still unresolved problem of how to construct 
the LSDA+D functional in order to avoid the double-counting problem 29]. Taking into 
account the above controversy, it is crucially important to treat the electron correlations (if 
any) in the most unambiguous manner. In the present work, we will try to pursue this strat¬ 
egy, first, by constructing the realistic model for Cr 02 and deriving all the parameters from 
first-principles calculations and, then, by solving this model within the dynamical mean-held 
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theory (DMFT), supplemented with the exact diagonalization (ED) method for the quan¬ 
tum impurity problem. We will show that stability of the HM FM ground state in CrC >2 is 
a highly nontrivial problem. If one considers only static electron correlations in the frame¬ 
works of either unrestricted Hartree-Fock (HF) or static DMFT techniques, the FM ground 
state can be formally obtained already in the minimal model, including only the closest to 
the Fermi level t- 2 g bands. However, the dynamic correlations tend to destabilize this state. 
Therefore, in order to explain the experimentally observed ferromagnetism in CrC> 2 , it is 
crucially important to consider other magnetic interactions and we argue that these are the 
direct exchange interactions between Wannier functions centered at different Cr sites and 
the polarization of the oxygen 2 p band. Another static approach - the so-called optimized 
effective potential (OEP) method, treating the correlation interactions in the random-phase 
approximation, produces a curious but unphysical insulating solution and further suppresses 
the tendencies towards ferromagnetism. This again emphasizes the importance of consistent 
treatment of correlation interactions in CrC> 2 . 

The rest of the article is organized as follows. In Sec. [H] we will explain details of our 
method: the construction of effective low-energy model (Sec. Hi]), the solution of DMFT equa¬ 
tions ISec. Ill Bl) . and the difference between unrestricted HF and static DMFT techniques. 
In Sec. IIIII we will present our results for interatomic exchange interactions and discuss 
them in many details: the (AE)^ 1 expansion for nearest-neighbor (NN) and next-NN in¬ 
teractions (Sec. IIII AD . the magnetic state dependence of interatomic exchange interactions 
(Sec. IIII Bl) . the behavior of long-range interactions (Sec. IIII CD . the contribution of the di¬ 
rect exchange interactions and the oxygen states (Sec. IIII DP . as well as results of the OEP 
method (Sec. IIII El) . Finally, in Sec. IIV1 we will present a summary of our work. 


II. METHOD 


A. Parameters of effective low-energy model 


In this section, we briefly remind the reader the main ideas behind the construction of 
effective low-energy model and present results of such construction for Cr02- The method¬ 
ological details can be found in the review article 30]. All calculations have been performed 
using parameters of experimental rutile structure (the space group P4 2 /rnrirri = D\fy 31]. 
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The model Hamiltonian, 



ij a ab i aa' abed 


(i) 


is formulated in the basis of Wannier orbitals {0 ia }, which are constructed for the mag¬ 
netically active Cr t- 2 g bands near the Fermi level, starting from the band structure in the 
local-density approximation (LSDA) (Fig. [T]). Here, cr(cr') = f /l are the spin indices, while 
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FIG. 1. (Color online) Total and partial densities of states of CrC >2 in the local density approxi¬ 
mation. The shaded light (blue) area shows the contribution of the Cr 3d states. The positions of 
the main bands are indicated by symbols. The Fermi level is at zero energy (shown by dot-dashed 
line). 


a, b, c, and d label three t 2 g orbitals, which have the following form in the global coordi¬ 
nate frame: |1) = ±^\xy)+ 2 Y-\3z 2 — r 2 ), \2) = -^\yz)±\zx), and |3) = \x 2 —y 2 ), where the 
upper and lower signs stand for the Cr site 1 and 2, respectively (see Fig. [2]). These orbitals 
are sometimes denoted as, respectively, \xy), \yz—zx), and \yz+zx), referring to the local 
coordinate frame [l2]. It is important that at Cr sites all three orbitals belong to different ir¬ 
reducible representations of the point group mmm = Fbh, meaning that all local quantities, 
including the crystal field, DMFT self-energy, and local Green’s function, will be diagonal 
with respect to these orbital indices. Moreover, the diagonal matrix elements will be the 
same for the Cr-sites 1 and 2. 
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FIG. 2. (Color online) Atomic electron densities, explaining relative positions of Cr t 2 g orbitals at 
the sites 1, 1', and 2. The oxygen atoms are indicated by the green spheres. 


Each lattice point % (j) is specified by the position r (r') of the Cr site in the primitive 
cell and the lattice translation R. Hence, the basis orbital (p ia ( r) = 0 Ta (r—R— t) is centered 
in the lattice point (R+r) and labeled by the indices r and a. The Wannier functions were 
calculated using the projector-operator technique 30, [32J and orthonormal linear muffin-tin 
orbitals (LMTO’s) 33H35I] as the trial wave functions. In physical terms, LMTO can be 
viewed as the localized atomic-like Wannier function constructed for the whole region of va- 
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lence states. Therefore, the projector-operator technique allows us to generate well localized 
Wannier functions for the t 2 g bands, that is guaranteed by the good localization of LMTO’s 
themselves. Then, the one-electron part of the model is identified with the matrix elements of 
LDA Hamiltonian (Hlda) in the Wannier basis: f“ fc T , +R = (0 Ta (r—r)|'HLDA|0r , fe( r— 

Since the Wannier basis is complete in the low-energy part of the spectrum, the construction 
is exact in the sense that the band structure, obtained from f“ 6 T , +R , exactly coincides with 
the one of LDA. 

The site-diagonal part of t l3 = [t“ 6 ] describes the crystal field splitting. It has the following 
form (in meV): 


/ 


\ 


/ 


-246 0 0 

hi = 0 60 0 

v 0 0 186 

The matrices of transfer integrals in the bonds 1-1' and 1-2 are given by 

—67 0 0 ^ 

tw = 


( 2 ) 


V 


-67 0 0 

0 -191 0 

0 0 158 


( 3 ) 




and 


ti 2 — 


/ 


V 


-15 0 
-28 0 
0 194 


0 

0 

-119 


\ 
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( 4 ) 


respectively. Other transfer integrals are considerably weaker [36]. The obtained values are 
in reasonable agreement with results of previous calculations jl2[. One interesting aspect is 
the large matrix element tff, = 158 meV, which is formally of the ddd type (see Fig. [21) and, 


therefore, supposed to be weak 


37]. Nevertheless, such large transfer integrals are possible 


due to pecu 


iar geometry of the Cr 02 lattice and contributions of intermediate O 2 p and Cr 


e g states [l2]. In terms of the Wannier functions, this means that the functions should have a 
sizable tail spreading to the oxygen and other Cr sites [l2]. Thus, already from this fact one 
can expect appreciable direct exchange interactions, which will be evaluated in Sec. IIII D1 
Another interesting aspect is the large asymmetric contribution t 3 \ = 194 meV, caused by 
the electron transfer via intermediate oxygen atom (see Fig. [2}i). The same mechanism is 
responsible for finite However, it is considerably smaller than tf^- 
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Matrix elements of the on-site Coulomb interactions can be also obtained in the Wannier 


basis as 


U ab cd= I dr J 

where the screened interaction u scr (r, r') is computed in the constrained random-phase ap¬ 
proximation (RPA) 38]. Since RPA is very time consuming, we apply additional approx¬ 


imations, which were discussed in 


30]. Namely, first we evaluate the screened Coulomb 


and exchange interactions between atomic Cr 3 d orbitals, using fast and more suitable for 
these purposes constrained LDA technique. After that, we consider additional channel of 
screening caused by the 3d —> 3d transitions in the polarization function of constrained RPA 
and project this function onto the 3d orbitals. The so obtained parameters of Coulomb in¬ 
teractions are well consistent with results of full-scale constrained RPA calculations without 
additional approximations. 

The obtained matrices of the on-site Coulomb interactions were fitted in terms of two 
Kanamori parameters 39]: the parameter of intra-orbital Coulomb interaction U = 2.84 eV 
and the exchange interaction J = 0.70 eV. The third Kanamori parameter - the so-called 
inter-orbital Coulomb interaction - can be obtained from U and J as IA' — U — 2 J. 


B. Dynamical mean-field theory 


Solution of the low-energy model, represented by Hamiltonian (JTj) , is a complicated nu¬ 
merical and methodological problem. In general, microscopic properties of a periodic mag¬ 
netically collinear system can be expressed via one-electron Green’s function k), 

which, in the reciprocal space, can be formally related to the frequency- and momentum- 
dependent self-energy £74 , k): 


G'T’^k) 


to — t(k) — ££4(0^ k) , 


(5) 


where t( k) is the one-electron part of the Hamiltonian (0Q) in the reciprocal space and all 
kind of correlation effects are described by £74^ k). 

The basic approximation, underlying the dynamical mean-held theory (DMFT), is that 
the self-energy is assumed to be independent on k: 


£7’V,k) ~S W H, 


( 6 ) 
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which becomes exact in the limit of infinite dimensions (or coordination numbers) 40]. The 
main idea of DMFT is to map the initial many-body problem for the crystalline lattice 
onto the quantum impurity one, surrounded by an effective electronic bath, and find self- 
consistently the parameters of this bath. Namely, the local (or site-diagonal) Green function 
of the crystal is given by 

G'T’V) = J]G U (w,k). 

k 

It can be further used to obtain the bath Green function, Q{l j), from the Dyson equation: 

(7) 


Then, new (A^u;) is obtained by solving the Anderson impurity model. The correspond¬ 
ing Hamiltonian is given by 


liimp ^2 E a.dl a daa + „ ^ Uabcddld ca , didder' 


abed, cr, a' 


'y 1 \Vapd\ a c p(T + H.c. + y ^ e 


pC\> a Cp<Tl 


ap,cr 


p,a 


( 8 ) 


where d(cr) and c(c^) are the electron annihilation(creation) operators for the impurity and 
bath states, respectively, V ap is impurity-bath hybridization, and E a ( e p ) are the noninter¬ 
acting energy levels of the impurity (bath). In order to obtain parameters of the Anderson 
impurity model, we adapt the following analytical form of the bath Green function (sepa¬ 
rately for each f 2s orbital a): 


= (u-E a -J2 — 

\ - (JJ — 

\ P 

and fit it in terms of E ai e p , and V ap . Generally speaking, the number of bath states p is 
infinite. However, in order to handle this problem numerically by means of ED, we discretize 
the bath and use a Enite number of electronic orbitals N s . It enables us to numerically 
diagonalize the impurity Hamiltonian (J8]) and obtain which is further identified with 
Then, using G a {u) = ^ 3 ( w ), the new self-energy can be found from the Dyson 
equation ©• After that, it is substituted into Eq. ([5]) to obtain new CH-t^k), and the 
problem is solved self-consistently. 

The ED method allows us to hnd the ground state as well as the low-lying excitations of 
the quantum impurity model. The standard numerical algorithms to treat the eigenproblem 
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are based on the matrix-vector multiplication, where the initial vector, matrix and net 
vector are kept in computer’s random-access memory (RAM). In order to make our model 
treatment realistic, it is necessary to take total N s ranging from 15 to 18. However, it would 
lead to the Hamiltonian matrix (jS]) of the dimensionality rsj (10 10 xl0 10 ), which makes the 
diagonalization procedure troublesome, even for modern multiprocessor computers. 

In this study we use the newly developed numerical ED scheme, based on the standard 
Arnoldi algorithm implemented in ARPACK program package 41], where the Hamiltonian 
matrix is not stored in the RAM, but efficiently recalculated ‘on the fly’ at each matrix- 
vector multiplication step. It makes the computational time to increase by only 10-15%. 
However, the amount of necessary RAM is decreased by 80-90%, giving us the possibility 
to perform realistic calculations with the large number of effective orbitals N s . Particularly, 
using this scheme, we were able to include 4 bath states per each Og orbital in the framework 
of DMFT. We have confirmed that the obtained electronic structure is well converged de¬ 
pending on the number of the bath states. The numerical calculations have been performed 
for the temperature T = 232 K, which is substantially larger than the magnetic transition 
temperature, and in the external magnetic field /j>bH = 5 meV, which is required in order 


to lift the magnetic degeneracy of multiplet states 


42], The example of electronic spectrum 


is shown in Fig. [3j which is in remarkable agreement with results of the previous DMFT 


studies 


13]- 


C. Static DMFT versus unrestricted Hartree-Fock approach 


The asymptotic high-frequency behavior of XflRu;) in DMFT is given by 43]: 


Ej(oo) — (U — 3 J)(nl + nl) +Un\ + (U - 2J)(n\ + nj), 


(9) 


where {rRl} are the self-consistent populations in DMFT: 


= -Im 

7T 


f£ F 


dwGj'V). 


Other matrix elements of XffiRoo) can be obtained from Eq. (J9]) by permutation of the spin 
and orbital indices. E^Roo) has the same form as the potential matrix in the unrestricted 


Hartree-Fock method 


30], but with different populations: in DMFT, these populations 


include the effect of frequency-dependence of the self-energy, while in HF, they do not. For 
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HF DMFT 



FIG. 3. (Color online) Partial densities of states as obtained in the unrestricted Hartree-Fock 
approach (left) and the dynamical mean-field theory (right) for the ferromagnetic state. The 
Fermi level is at zero energy (shown by dot-dashed line). 


TABLE I. Self-consistent orbital populations {nl’'*'} for the ferromagnetic states, as obtained in 
DMFT and unrestricted Hartree-Fock calculations for the low-energy model. 



t 

n[ 

t 

n 2 

t 

n 3 

1 

n\ 

4 

n 2 

4 

DMFT 

0.934 

0.587 

0.431 

0.010 

0.020 

0.023 

HF 

0.999 

0.710 

0.291 

0 

0 

0 


the FM state, these populations are summarized in Table [H Besides small population of the 
J,-spin states (and, therefore, small deviation from the HM behavior), the main difference 
between DMFT and HF is in the orbital polarization of the "t-spin states (see Fig. [3]). The 
first orbital is practically fully occupied in both approaches. The population of other two 
orbitals tends to be nearly equal in DMFT, while in HF these states are strongly polarized 
and there is an additional redistribution of electrons between v} 2 to rig. Finite values of 
is the natural result of DMFT calculations for the HM magnets^ which is related to 
the existence of nonquasiparticle J,-spin states near the Fermi level [7, 13]. Nevertheless, we 
have found that these states have a minor effect on the behavior of interatomic exchange 
interactions and, from this point of view, CrC >2 can be treated as HM ferromagnet, even in 
DMFT. 
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III. INTERATOMIC EXCHANGE INTERACTIONS 


We consider the mapping of electron model (JTJ) onto Heisenberg model with S — 1: 




j+i- 


In these notations, J* is the exchange coupling between two Cr sites, located in the origin 
(0) and in the point i of the lattice, relative to the origin. The mapping onto the spin 
model implies the adiabatic motion of spins when all instantaneous changes of the electronic 
structure adjust slow rotations of the spin magnetic moments. The parameters of this model 
can be obtained by using the theory of infinitesimal spin rotations 4dJ, l45i] : 


1 


f£F 


Ji = ^Im / doj Tr L \ A£M(%(w)AE(w)Gk(a;) 


( 10 ) 


where Gg^ca) = [cu — t — i s the one-electron Green function between sites 0 and 

i, AI = 5T — id and Tr j j denotes the trace over orbital indices. The parameters {Jj} given 
by Eq. (gnD are nothing but the second derivatives of the total energy with respect to the 
rotations of spins. Therefore, this definition of the Heisenberg model is valid only for small 
rotations of the magnetic moments near the FM state and characterizes the local stability 
of this state. The effect of finite rotations will be discussed in Sec. mm 

The parameters of interatomic magnetic interactions, obtained in the theory of infinitesi¬ 
mal spin rotations, are listed in Table UTI and their behavior is explained in Fig.Q] We note 
the following: (i) As expected, the FM ground state is stabilized by the NN and next-NN in¬ 
teractions (Ji and J 2 , respectively). The values of these interactions, obtained in DMFT and 
unrestricted HF approach, are surprisingly close, while static DMFT overestimates both of 
them; (ii) Besides strong FM interactions J\ and J 2 , there are several types of AFM interac¬ 
tions, operating in the 5th, 6th, 7th, and 8th coordination spheres, which tend to destabilize 
the FM state. These interactions are especially strong in the case of DMFT. In the next 
sections, we will elucidate the microscopic origin of such behavior and its consequences on 
the properties of Cr0 2 . 


A. Double exchange and beyond 


The ferromagnetism of Cr0 2 is frequently attributed to the DE mechanism 
This is the very important point, which needs to be clarified. 
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TABLE II. Parameters of interatomic exchange interactions (in meV) as obtained in the theory 
of infinitesimal spin rotations with three different types of approximations for the self-energy: the 
unrestricted Hartree-Fock approximation (HF), the dynamical mean-field theory (DMFT), and 
the static limit for the DMFT self-energy E(u;—»oo) (SDMFT). The notations of parameters Jj are 
explained in Fig. [4j 


parameter 

HF 

DMFT 

SDMFT 

Ji 

14.06 

16.35 

19.63 

h 

12.26 

12.14 

13.65 

J 3 

1.16 

0.60 

0.82 

■h 

0.96 

0.35 

0.78 

Jb 

-0.39 

-1.15 

-0.72 


-1.87 

-1.85 

-1.78 

J< 

-1.21 

-2.58 

-1.45 

J> 

-3.26 

-4.19 

-3.25 


-0.31 

-0.94 

-0.47 


-0.46 

-2.44 

-1.00 


In the HM regime, all poles of Gh(u;) are located in the unoccupied part of the spectrum 
and AS can be regarded as a large parameter. Note that the existence of the small weight 
of nonquasiparticle j,-spin states near the Fermi level does not alter this conclusion, which 
will remain true even in the case of DMFT. This justifies the use of the (AS) -1 expansion 
in the occupied part 22]: 


GV) = -AE-‘Xi([AEGf]- 1 )", 

n =0 

which follows from the identity — AS^ . The n = 0 term of this expansion 

contains only site-diagonal elements and, consequently, does not contribute to Eq. m- 
Therefore, in the HM state, can be presented as an infinite series: 


■>i = E 47 (H) 

71=1 
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FIG. 4. (Color online) (Left) Distance dependence of interatomic exchange interactions as obtained 
in the theory of infinitesimal spin rotations with three different types of approximations for the 
self energy: the unrestricted Hartree-Fock approximation (HF), the dynamical mean-field theory 
(DMFT), and the static limit for the DMFT self-energy E(cj— >oo) (SDMFT). (Right) Lattice of 
Cr sites with the notation of interatomic exchange interactions. 


where 


4 ” ) = “Im r dwit L {g;,m ([AE(uOGV)]- 1 )” aeh} . (12) 


The n — 1 term corresponds to the DE interaction, which can be easily found analytically 

22 ]: 

1 C F r . -1 

■io \ • (13) 


4 1} = -^ Im / duTr L <( G^( u)u 


Moreover, using the identity G'(uj) cu — t — £l(u;) 


= 1, it is straightforward to show that 


E-4 11 = -1^, 


(14) 


where E^ n is the kinetic energy (per one Cr site): 


E kin = --Im [ du Tr L {G'J 0 (w) u - t 00 - Xf^a;) 
71" .1 I 


From Eq. (ITT) it is clear that the main interactions J- 1 ^ should be positive (or ferro¬ 
magnetic). This equation is nearly perfectly reproduced by our calculations. For instance, 
we have obtained the following values of DE interactions in DMFT: J ^ = 28.95 meV, 

J 2 — 19.44 meV, = 1.33 meV, and = 1.59 meV. Since the Wannier functions are 
localized and the transfer integrals connecting more remote sites are small, the corresponding 
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to them parameters of DE interactions are also small. Then, by considering the sum of DE 
interactions up to the fourth coordination sphere, we find 2= 231.50 
meV, which readily reproduces 99 % of the kinetic energy — \E^ n = 234.13 meV. A very 
similar conclusion holds in unrestricted HF and SDMFT calculations. 

Moreover, SDMFT yields very similar parameters of the main DE interactions: jj 1 = 
28.71 meV and = 19.82 meV, which are practically undistinguishable from the ones in 
DMFT. However, the parameters obtained in the HF calculations are considerably smaller, 
especially for the nearest neighbors: j [ 1) = 23.75 meV and = 19.03 meV. Such behavior 
is directly related to the orbital polarization and additional splitting of the states 2 and 3 
around the Fermi level (see Fig. |3j which tend to decrease \E^ in \ and, therefore, the values 
of DE interactions. For instance, in the case of jj 1 ^, all transfer integrals fl3]) are diagonal 
with respect to the orbital indices. Therefore, as the orbitals 2 and 3 become, respectively, 
more and less populated in the case of HF calculations (see Table [I]) , the DE interaction 
j \ 1 will decrease. In the case of the situation is less straightforward, because the 
transfer integrals (J4J) mix the orbitals 2 and 3, counterbalancing the change of the orbital 
occupations. 

Other contributions to can be found numerically. Particularly, Jj 2 ' 1 is of the Erst order 
of [AS(w)] _1 . It contains the contributions of superexchange interactions and the exchange 
processes between sites separated by two hoppings. In static case, all these parameters can 
be expressed via moments of the local density of states [21], However, in dynamic case such 
simple relationship does not exist, because of the frequency-dependence of E. 

The results of these calculations are shown in Fig. [5] Both J { 1 " ) and display some 
characteristic oscillating behavior, where the odd FM contributions are partially compen¬ 
sated by the even AFM ones. This tendency is observed in all the calculations, based on the 
unrestricted HF, SDMFT, and DMFT techniques. The main difference is the convergence 
of Y, n J i n \ which is noticeably slower in DMFT: the frequency-dependence substantially 
reduces Re[AE] in the occupied part, especially in the region close to the Fermi level (see 
Fig. |6]) and, therefore, slows down the convergence of the (AE) -1 expansion. On the other 
hand, Im[AE] is relatively small in the occupied part and does not play a significant role. 
Another important aspect is the cancelation of FM and AFM contributions to J\ and J 2 - 
As was discussed above, the unrestricted HF approach yields somewhat weaker FM DE 
contributions ./[ 1 and ./(j , due to the orbital polarization effects. However, the next AMF 
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FIG. 5. (Color online) Results of the (AS) -1 expansion for the nearest-neighbor (top) and next- 
nearest-neighbor (bottom) exchange interactions. The individual contributions, j\ n \ are shown on 
the left panel and their sum - on the right panel. The asymptotic values of J\ and J 2 are shown 
by the dash-dotted lines. 


contributions and ' are also weaker due to the larger spin splitting AS in comparison 
with DMFT. Thus, the total values of J\ and J 2 , obtained after summation of all these 
contributions, appear to be very close in the case of HF and DMFT. 


The series J^ n> is practically converged for n — 5, where these sums are close to the 
saturated values of -J\ and J 2 . The major FM contribution to J\ and J 2 is indeed due to the 
DE mechanism (n = 1). However, this contribution is not the only one and, at least, the 
n = 3 term is also very important in stabilizing the FM ground state. Thus, already from 
this point of view, it is not quite right to consider Cr0 2 as the DE system: the behavior of J\ 
and J 2 involves other important mechanisms besides the double exchange and superexchange 
interactions, which are considered in the conventional DE model If], 20], 
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FIG. 6. (Color online) Frequency dependence of intraatomic spin splitting AS(w) in DMFT. The 
static limit AS(w->oo) is shown by dashed lines. The Fermi level is at zero energy (shown by 
dot-dashed line). 

B. Magnetic-state dependence of interatomic exchange interactions and Curie 
temperature 


Generally speaking, the exchange interactions (TTOT) depend on the magnetic state in which 
they are calculated. This dependence reflects the change of the electronic structure in differ¬ 
ent magnetic states and such information is incorporated in the one-electron Green function 


GffA(u;). The magnetic state-dependence of exchange interactions may 
ical origin. For instance, it can be the orbital ordering in insulating 


rave different phys- 


46] or metallic 47] 


systems, or simply the change of the bandwidth in metallic compounds depending on the 
magnetic state 19]. The theory of infinitesimal spin rotations 4d], [45] is more suitable for 
the description of effects, which are related to small variations of magnetic moments near 
the ground state (for instance, the spin waves). Generally speaking, it is not applicable for 
the analysis of large perturbations, such as the spin disorder near the Curie temperature 
(T c ), unless the exchange interactions do not depend on the magnetic state. 
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The comparison of interatomic exchange interactions, calculated in the FM and AFM 
states using the theory of infinitesimal spin rotations in the case of HF and DMFT tech¬ 
niques, is given in Table IIIII (throughout this work we consider the simplest AFM configu¬ 
ration, where the corner and body-centered Cr moments in the single unit cell are oriented 
antiferromagnetically). One can clearly see that the exchange interactions are quite sensitive 

TABLE III. Parameters of interatomic exchange interactions (in meV) and corresponding Curie 
temperature (in Kelvins), obtained using different techniques and starting conditions, such as the 
theory of infinitesimal spin rotations near the ferromagnetic (F) and antiferromagnetic (A) state in 
the frameworks of unrestricted HF and DMFT methods, as well as the mapping of the total energies 
obtained in the self-consistent HF calculations for the spin-spiral configurations onto Heisenberg 
model (SCHF). The notations of parameters Jj are explained in Fig. [4j The dash sign in the row 
Tq means that for the given set of parameters the ferromagnetic state is unstable. 


parameter 

HF 

DMFT 

SCHF 

F 

A 

F 

A 

Ji 

14.06 

23.77 

16.35 

18.14 

11.00 

J 2 

12.26 

14.91 

12.14 

6.73 

15.43 

J 3 

1.16 

0.25 

0.60 

0.21 

2.90 

Ja 

0.96 

1.39 

0.35 

-1.08 

1.43 

J 5 

-0.39 

1.39 

-1.15 

-2.66 

0.10 

•h 

-1.87 

-0.14 

-1.85 

0.22 

-1.36 

J< 

-1.21 

-6.21 

-2.58 

-2.68 

-4.13 

J> 

-3.26 

-7.99 

-4.19 

-5.03 

-4.13 


-0.31 

-0.94 

-0.94 

-3.57 

-1.41 

Ji 

-0.46 

-3.05 

-2.44 

-1.60 

-1.41 

T C 

581 

— 

— 

— 

684 


to the magnetic state in which they are calculated. Generally, the AFM structure remains 
unstable and is not the ground state of CrOo. Nevertheless, the AFM spin alignment tends 
to reconstruct the electronic structure (Fig. [7]) in such a way as to additionally stabilize the 
FM interactions J\ in the NN ferromagnetically coupled bond. Moreover, in DMFT, the 
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FM interactions J 2 in the antiferromagnetically coupled next-NN bond is strongly reduced, 
that also works in the direction of stabilizing the AFM state (and destabilizing the FM one). 
Similar tendency holds for longer-range interactions. Thus, if one tries to use the parameters 



FIG. 7. (Color online) Partial densities of states as obtained in the unrestricted Hartree-Fock 
approach (left) and the dynamical mean-field theory (right) for the antiferromagnetic state. The 
Fermi level is at zero energy (shown by dot-dashed line). 


obtained in the AFM state in order to describe the FM state, one can easily find that this 
FM state will be unstable, even in the unrestricted HF approach. Perhaps, this was an 
extreme example, and below we will consider a more realistic strategy for the evaluation of 

T c - 


Taking into consideration the strong magnetic state-dependence of exchange interactions, 
we tried to go beyond the theory of infinitesimal spin rotations and evaluated the exchange 
interactions using results of self-consistent total energy calculations for spin-spiral config¬ 
urations with arbitrary wavevectors q. Namely, using generalized Bloch theorem 48], we 
performed the unrestricted HF calculations for spin-spiral configurations, where the direc¬ 
tions of magnetic moments varied as 


e r+R 


^ cos(t + R) • q ^ 
sin(r + R) ■ q , 


0 / 


calculated the total energy (A q ) for each q, and evaluated the exchange interactions as the 
Fourier transform of E n . The results are also listed in Table m in the column ‘SCHF’. 
Particularly, we expected that the exchange interactions, obtained by mapping the total 
energies of the spin-spiral configurations onto the Heisenberg model, should provide a good 
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estimate for Tq. The latter was evaluated using Tyablikov’s random-phase approximation 


49]. The results are also listed in Table Hill and can be summarized as follows: as long as we 


use the unrestricted HF approximation, supplemented either with the theory of infinitesimal 
spin rotations near the FM state or with the self-consistent spin-spiral calculations for finite 
q’s, running through the first Brillouin zone, Tq is even overestimated in comparison with 
the experimental data, meaning that the FM state is indeed very robust. However, when we 
switch to more rigorous DMFT technique, the FM state appears to be unstable because of 
the longer-range AFM interactions (and any numerical estimates of Tq in this case become 
meaningless). This is a very serious problem, which we will discuss in details in the next 
section. 


C. Long-range interactions and stability of the ferromagnetic state 


In Sec. 1III A[ we have seen that, as far as the NN and next-NN interactions are con¬ 
cerned, unrestricted HF and DMFT techniques very similar results. Nevertheless, there is 
an important difference in the behavior of longer-range interactions, which has fundamental 
consequences. Since the frequency-dependence reduces intraatomic spin splitting Re[AE(o;)] 
near the Fermi level, the series (fTT|) converges somewhat slower in the case of DMFT. Besides 
oscillating behavior depicted in Fig. [5j smaller Re[AE(u;)] is responsible for larger spatial 
extension of the exchange interactions. This can be directly seen from the construction 
^[AE(o;)G 1 '(a;)]~ 1 j in Eq. (fT2j) : since [^(o;)]^- 1 = t t j for and the transfer integrals 
are typically restricted by only few coordination spheres, the n-order term will include the 
processes, which connect two remote sites 0 and i by n sequential hoppings between near¬ 
est or next nearest neighbors. Obviously, such contributions will be stronger for smaller 
Re[AE(u;)]. Moreover, the number of nodes of the integrand in Eq. (fT21) increases with the 


distance between 0 and i 


50 


51]. Therefore, it is possible that some of these long-range 


interactions can easily become antiferromagnetic. Such behavior is clearly seen in Table HH 
and Fig. [4j besides FM interactions, there are several relatively strong AFM interactions, 
connecting the sites in the 5-8 coordination spheres. These interactions are stronger in 
DMFT because of smaller spin splitting Re[AE(u;)]. 

The appearance of AFM interactions naturally rises the question about stability of the 
FM state and whether it is indeed the magnetic ground state of the considered model. In 
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order to investigate this problem, we evaluate the spin-wave dispersion, u;(q), using the 
interatomic exchange interactions obtained in the theory of infinitesimal spin rotations. In 
the PA-z/mnm structure, containing two magnetic sublattices, u;(q) can be obtained from 
the diagonalization of the 2x2 matrix (for S'=l): 


fi(q) 


Jn(q) — Jo ^12 (q) 
^21 (q) ^22 (q) — Jo 


) 


where J a p{ q) is the Fourier image of magnetic interactions between sublattices a and /3, and 
J 0 = Jn(0) + J 12 (0) . In principle, due to the symmetry properties, J 22 (q) can be related to 
Jn(q*) in some other q-point. The same holds for Ji2(q) and J 2 i(q*)- The results of these 
calculations are shown in Fig. [HI The negative spin-wave frequencies signal that the FM state 
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FIG. 8. (Color online) Results of calculations of the spin-wave dispersion with the parameters, 
obtained in the theory of infinitesimal spin rotations in the case of HF, SDMFT, and DMFT 
techniques. Notations of the high-symmetry points of the Brillouin zone are taken from Q. 

is unstable. One can clear see that as long as we use the static HF and SDMFT techniques, 
there is no problem with the stability of the FM state, and Tq is even overestimated in 
comparison with the experimental data (Table IIIID . Thus, one could naively think that the 
FM state is very robust. Nevertheless, in DMFT, which is definitely the most rigorous 
approach among the considered ones, the FM state appears to be unstable. This instability 
occurs along three high symmetry directions of the Brillouin zone (T-X, T-M, and T-Z). This 
is a very serious problem, meaning that there should be additional factors, which are not 
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taken into account in the low-energy electron model for the t 2g bands and which stabilize 
the FM state. This problem will be studied in the next section. 


D. Direct exchange interactions and contributions of the oxygen states 


In this section, we evaluate the change of magnetic energy, caused by the polarization of 
the O 2 p band and other contributions, which are not taken into account in the minimal 
model for the t 2g bands. 

For these purposes, after solution of the low-energy model in DMFT, we go back from 
the Wannier basis {< f > Ta } of the model to the original LMTO basis {Xvb}- 

0ra(r - r) = VraXvbir - v), (15) 


vb 


and construct the spin magnetization density, m(r) = ny(r)—n^(r), associated with the Cr 
t 2g band. This m(r) has major contributions at the Cr sites as well as some hybridization- 
induced contribution at the oxygen sites. Following the philosophy of the low-energy model 
j3o| . the interaction of m( r) with the rest of the electronic states should be well described 
already at the LSDA level. Therefore, our strategy is to evaluate, in LSDA, the exchange- 
correlation (xc) held b( r) = ry(r) — ty(r) (tyy being the xc potential in LSDA), which is 
induced by m( r) and polarizes the O 2 p band, and find the self-consistent change of m( r) 


and b(r), caused by the interaction between t 2g and O 2 p bands. For 


convenient to use the self-consistent linear response (SCLR) theory 53]. For simplicity, let 


hese purposes, it is 


us consider the discrete lattice model and assume that all weights of m( r) are concentrated 
in the lattice points: m( r) = m v S(r — v), where m v is the local magnetic moment at the 
site v. Furthermore, we recall that LSDA is conceptually close to the Stoner model, where 


the xc energy is given by 


54j]: 


E xc — ^ ^ I v m v . 


(16) 


In practical calculations, the parameters {I v } can be found using the values of intraatomic 
spin splitting and local magnetic moments in LSDA. Meanwhile the intraatomic spin split¬ 
ting itself can be obtained using LMTO parameters of the centers of gravity for the 'f- and 
J,-spin states 34], which yields Iqi = 0.98 eV and Jo = 1.68 eV. 

Then, the self-consistent held can be found as 


b = 


1 +11Z 


-i 


5°, 
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where we have introduced the vector b = [&„] and the tensors X = [I V 5 VV >] and 71 = [7 Z vv >]. 
In this equation, b° = Xrh is the xc field induced by the Cr t 2g band, and the response tensor 
7 Z is obtained in the first order perturbation theory for the wavefunctions, starting from the 
nonmagnetic LDA band structure: 


occ unocc BZ 


ttw-VVVV 


(s~iva\*/^va (s~yv'b\*sxv'b 
l°nkJ ^n'kK^n'k) °nk 


+ c. c. 


(17) 

7 , , l ^nk &n f k J 

ab n n' k v 

where {C^g} are the coefficients of the expansion of the LDA wavefunctions over LMTO’s, 
{e n k} are the LDA eigenvalues, and k runs over the first Brillouin zone (BZ). Moreover, 


similar to the constrained RPA 


38], we have to exclude from Eq. (fTTj) contributions, where 


both indexes n and n' belong to the Cr t 2g band. In the perturbation theory, such terms 
describe the change of the magnetization in the t 2g band, which are caused by the LSDA 
potential. However, in the low-energy model, this part is replaced by the more rigorous 
DMFT solution with the screened Coulomb interactions. Therefore, in order avoid the 
double counting, such contributions should be excluded in the process of SCLR calculations. 
In practice, n runs over the occupied O 2 p bands and n' runs over the unoccupied Cr t 2g 
and e g bands. 

Once the self-consistent field b is known, the change of m and b, caused by the polarization 
of the oxygen band, can be found as Sm = —7 Zb and 5b = 15m, respectively. Since the 
O 2 p band is occupied, the net change of magnetic moment will vanish: Y2 V $ m v — 0, 
irrespectively on the type of the magnetic order. Nevertheless, the individual moments 5m v 
can be finite and contribute to the total energy. The corresponding energy change, caused by 
the magnetic polarization of the oxygen band, consists of two parts: 5E po1 = 5E^f_ 0 + 5E^ 1 , 
where 5E^f_ 0 = 5m T Xm is the interaction of 5rn v with the “external” xc field, created 
by the Cr t- 2g band, and h7?o o1 is the energy change caused by 5m in the O 2 p band. It also 
consists of two parts: 5E^ 1 = 5E sp + 5E dc , where 5E sp is the single-particle energy, which 
can be found in the second order of 5b as 5E sp = \5b T 7Zb 


53], and 5E dc = jbfh^Xbrh is the 


double-counting energy, where 5rh T is the row vector, corresponding to the column vector 
5m. In all these calculations, it is assumed that the magnetic energy of the t 2g band itself 
is described by DMFT. 

The polarization energy 5E po[ may have different values in the case of the FM and AFM 
alignment of spins and, thus, contributes to interatomic exchange interactions. Below we 
evaluate this effect in Cr02- The magnetic moments are listed in Table HVl and the energies 
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are in Table [V] 


TABLE IV. The values of local magnetic moments at the chromium and oxygen sites {rn v } as 
obtained in DMFT calculations for the isolated t 2 g band in the case of ferromagnetic (F) and 
antiferromagnetic (A) alignment of Cr spins, and the moments caused by the polarization 

of the O 2 p band. All values are in pb- 




F 



A 



m v 


5m v 

m v 


5m v 

Cr 

1.628 


0.594 

1.392 


0.584 

0 

0.134 


-0.297 

0.029 


-0.089 


TABLE V. The energy changes (in meV per one formula unit), caused by the magnetic polarization 
of the O 2 p band in the ferromagnetic (F) and antiferromagnetic (A) states: the interaction energy 
between Cr t 2 g and O 2 p bands (AEq^q), the magnetic energy in the O 2 p band (AEq 01 ), and 
the total energy ( 5E po1 = q + <5Fq 01 ). All values were derived using DMFT magnetization 

density for the t 2 g band. 


F A 


r)F po1 

0 rj Cr-O 

-449.27 

-434.67 

5F po1 

99.04 

66.42 

5E po1 

-350.24 

-368.24 


The spin moments m v are redistributed between Cr and oxygen sites. As expected for 
the FM state, the total moment is mo + 2mo = 1.9 pb, which is totally consistent with the 
value obtained in the Wannier basis (see Table []]). The small deviation from 2 pb is caused by 
nonquasiparticle 4,-spin states near the Fermi level. The moments m v and 5m v are parallel 
at the Cr sites and antiparallel at the oxygen sites. This tendency is consistent with results 
of first-principles calculations and can be deduced from the analysis of hybridization between 
Cr 3 d and O 2 p states [55|]. Therefore, the negative sign of 5Ec r -o is due to the contributions 
of the Cr sites, which are partly compensated by positive contributions of the oxygen sites. 
The absolute value of 5E Cr -o is larger in the FM state, mainly because m<> and 5mc T are 
larger. Thus, the Cr-0 interaction additionally stabilizes the FM state. The contribution of 
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the O 2 p band to the magnetic energy is positive. This is because the 0 2 p band itself does 
not favor the magnetism and any magnetic polarization of this band will increase the total 
energy. This also explains why SEq is smaller in the AFM state: the magnetic moments 
Srri v are smaller and, therefore, the magnetic perturbation of the O 2 p band is also smaller. 
In Cr02, the second effect (SEq) dominates and the polarization of the O 2 p slightly favors 
the AFM alignment. The corresponding energy difference between FM and AFM states, 
A E po1 = 8E pol (F)—8E pol (A), is about 18 meV per one formula unit. 

Another contribution to the magnetic energy is related to the direct interactions between 
Wannier functions in the t 2g bands 26), |27j, which are centered at different Cr sites. They 
are not taken into account in the low-energy model, because the latter treats only on-site 
Coulomb and exchange interactions. Nevertheless, these interactions can be evaluated in 
LSDA. First, let us evaluate the difference of LSDA xc energies between FM and AFM states 
in the t 2g band, A E xc = 8E XC (F)—5E XC (A), using the values of magnetic moments {m v } from 
Table UV] and Eq. (fl6l) for E xc . This yields A E xc = —199.65 meV per one formula unit, where 
the main contribution (about 93%) comes from the Cr sites. This energy difference favors 
the FM alignment. Then, we note that the xc interaction between Wannier orbitals centered 
at the same Cr site is already taken into account in the low-energy model in the framework 
of DMFT. Therefore, we should subtract this on-site “self-interaction” (SI) part from the 
LSDA xc energy difference. This can be done as follow. Using the spin magnetization matrix 


<•£ F 


M T = [AC] = 


-Im 


doj 


7T 


Gf» - G*» 


obtained in the Wannier basis at the Cr site r, and the expansion (fT5]i over LMTO’s, we 
evaluate magnetic moments, that are produced by Al r at the central and neighboring to it 
sites v: 


m v = Y. («“)* 

abc 

The difference between m v and m v is that fh v is the contribution of the single Cr site r to 
the magnetic moment at the site v, while rn v takes into account the contributions of all sites 
of the Cr lattice. Therefore, rli v at the central site r is substantially smaller than m v (■ fh v — 
1.169 and 1.099 /i B in the FM and AFM state, respectively). The total moment 'Yh v fh v 
in the FM state is only 1.543 fi b, which also substantially deviates from Y2 v m v = 1-9 
/r B - Then, we evaluate the SI energy, which is also given by (fT6lh but after replacing 
{rn v } by {rh v }. This yields the following energy difference between FM and AFM states: 
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AEsi = §Esi(F)—6Esi(A) = —39.19 meV per one formula unit. Therefore, by subtracting 
the SI term, we will additionally shift the energy balance in the favor of antiferromagnetism. 

Thus, by combining all the contributions, the total energy difference A E = A E po1 + 
AE XC — AEsi is about —142.46 per one formula unit. By mapping this total energy difference 
onto the Heisenberg model and assuming that it contributes only to the next-NN interac¬ 
tions, one can find the following correction to this interaction, arising from the polarization 
of the oxygen band and direct exchange interactions in the t- 2 g band: AJ 2 = —AE/8 = 17.81 
meV. The spin-wave dispersion, which takes into account the additional FM contribution 
A J 2 is plotted in Fig. [9] in comparison with results of regular DMFT calculations for the 
isolated t 2 g band. One can clear see that all cu(q) in this case become nonnegative and the 



FIG. 9. (Color online) Results of calculations of the spin-wave dispersion with DMFT parameters 
obtained for the isolated t 2 g band (solid line) and after taking into account the additional ferro¬ 
magnetic contribution AJ 2 = 17.81 meV, arising from magnetic polarization of the oxygen band 
and direct exchange interactions in the t 2 g band (dotted line). Notations of the high-symmetry 


points of the Brillouin zone are taken from 


52]. 


FM state is stable. Thus, the magnetic polarization of the oxygen band and direct exchange 
interactions in the t 2g band play a very important role in the stability of the FM state in 
Cr0 2 . 
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E. Optimized effective potential method and importance of dynamic correlations 


In this section, we discuss results of the optimized effective potential method (OEP) 56 
60], which we consider mainly for pedagogical purposes, in order to emphasize the importance 


of careful treatment of the corre 
Sham density functional theory 


ation effects. OEP is a numerical realization of the Kohn- 


61] , where 


1. the one-electron band structure is obtained from solution of Schrodinger equations 
with some effective static local potential v: 


(tk ~P v) |Cnk) £nk|Cnk)) 


(18) 


2. the obtained band structure is used to calculate the total energy 

E = Ek in + Ec + Ex + E corr , (19) 

consisting of kinetic (E^ n , which also includes the energy of crystal-field splitting), 
Coulomb (Ec), exchange (Ex), and correlation (E corr ) parts; 

3. the parameters of effective potential v are found numerically, so to minimize the total 
energy (USD. 


Thus, the OEP method provides some alternative possibility for the construction of static 
potential, which, in addition to the standard Coulomb and exchange contribution, includes 
the effect of correlation interactions and, in this sense, can be regarded as a step beyond the 
HF approximation. 

This aforementioned OEP procedure was implemented for the solution of the low-energy 
model for the i 2fl band 62], Here, we assume that the one-electron band structure is half- 
metallic and all minority-spin states are unoccupied. Therefore, we drop the spin indices, 
but keep in mind that both the potential and electronic structure are referred to the '('-spin 
states. Then, because of the symmetry, the potential matrix is diagonal v = [v a b8 TT >8 a b\ and 
does not depend on the indices r = 1 or 2 of the Cr-atoms in the primitive cell. Therefore, 
the effective potential has only two independent parameters (apart from the constant energy 
shift): A 2 _i = V 22 —V 11 and A 3 _ 2 = v 33 — r> 22 . Note also that the eigenvector |c n k) in Eq. (TT8|) 
is the row-vector of the form |c„k) = [ c nkJ- 









The correlation energy can be evaluated in RPA as 
1 


63 


64J: 


Ecmr 4 n S 


duj Tr |In l — P(ioj,q)U l — UP(iuj,q) + 2P(iui, q)f/|, (20) 


where U is the matrix [U a b c d\ of the on-site Coulomb interaction and P = [-P^d] * s th e 
polarization in the imaginary frequency: 


occ unocc 


P TT ' (ii ) n'l — \ \ ' \ ' ^( g ' rak g n'k+q) q-r* 6r ct'* dr' 

r abcd\ LUJ i ~ / y y y 2 _i_ ( c x 2 Vk+q C nk C nk C n'k+q- 

„ “T W ~ £ n'k+qj 


( 21 ) 


The matrix multiplication in Eq. (|2Q|) implies the summation over two intermediate orbital 
indices: ( UP) T a l' cd = Y,e} U abefPlj' cd and the u -integration has been performed using 10- 
points Gaussian quadrature method [641. 

By applying this OEP approach, we expected that the correlation effects, beyond the 
HF approximation, will reduce the orbital polarization and yield an improved description, 
at least for the majority-spin states and DE interactions. Moreover, we expected RPA to 
work reasonably well for metallic systems, such as Cr0 2 . Since the RPA total energy of 
HM systems does not depend on the position of unoccupied J,-spin states, we cannot easily 
determine in the framework of this method the spin-splitting AS and the parameters of 
exchange interactions, which depend on AS and G^{u). Nevertheless, at least we should be 
able to evaluate the DE interactions, which do not depend on AS. 

However, somewhat surprisingly, we have obtained very curious, but unphysical result: 
the correlation interactions, treated in RPA with the static effective potential, tend to ad¬ 
ditionally stabilize the orbital ordering and increase the orbital polarization, leading to the 
insulating solution, which is shown in Fig. [TQI The reason for such unphysical behavior is 
that, even when the HF energy ERf = A’kin+A-c+A’x reaches its minimum, E corr continues 
to decrease as a function of A 3 _ 2 (Fig. [IT]). Obviously, E corr decreases when the polarization 
decreases (note that P is the negative-defined matrix in the imaginary frequency). Then, 
there are two competing effects. On the one hand, the additional splitting of orbitals 2 
and 3 across the Fermi energy is expected to suppress the correlation interactions. This is 
rather general property of correlation energy, which follows from the perturbation theory 
analysis 65]. If there were no transfer integrals, connecting the orbitals 2 and 3, the ef¬ 
fect of A 3 _2 would be equivalent to the scissors operator and the behavior of E corr would 
be totally described by the above mentioned mechanism (which indeed dominates for large 
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FIG. 10. (Color online) Electronic band structure of Cr02, obtained in the OEP approach: (Left) 
Total and partial densities of states of three t 2 g orbitals and (Right) band dispersion along high- 
symmetry directions of the Brillouin zone (notations of the high-symmetry points are taken from 
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FIG. 11. (Color online) Results of energy minimization in the OEP method versus the splitting 
in the potential between atomic levels 3 and 2: the Hartree-Fock part of the energy Ehf = 
in + Eq + Ex, the correlation energy E corr , and the total energy E = Ehf + E corr . 

A 3 _ 2 ). Nevertheless, the strong hybridization between orbitals 2 and 3 [see Eq. fllD] may 
change this canonical behavior. First, we note that, in order to produce a large contribution 
to E corr , one should activate the channels involving the large Coulomb matrix elements U aacc 
(where a and c are 2 or 3). This can be done only if the polarization matrix has sizable 
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elements of the same P aac c type [see Eq. (120|) ]. Such matrix elements are indeed produced by 
the hybridization effects in the insulating state [see Eq. f[2T]) ]: because of the hybridization, 
the orbital 3 may have a substantial weight in the occupied part of the spectrum (so as the 
orbital 2 in the unoccupied part), yielding finite matrix elements P aacc . Thus, we believe 
that the decreasing of E corr in Fig. [TT| is a specific property of Cr0 2 and related to the 
strong hybridization between occupied and unoccupied orbitals in the orbitally polarized 
state. Nevertheless, such a behavior is, of course, unphysical and this example demonstrates 
again the importance of explicit consideration of dynamic correlations. 

In order to conclude this section, let us evaluate the consequences of the exaggerated 
orbital polarization and gap opening on the DE interactions. The kinetic energy (without 
the energy of the crystal held splitting), obtained in the OEP method, is only —168.57 
meV per formula unit and corresponding parameters of DE interactions can be estimated as 
= 3.80 meV and = 8.74 meV. Thus, as expected, the DE interactions are strongly 
underestimated in the OEP approach. 

IV. SUMMARY AND CONCLUSIONS 

We have presented detailed theoretical analysis of interatomic exchange interaction in 
Cr0 2 , which was based on realistic low-energy model, derived from the first-principles elec¬ 
tronic structure calculations, and have involved various techniques for treating electron cor¬ 
relations in the narrow t 2g band, ranging from the static Hartree-Fock approximation to the 
dynamical mean-held theory. Such analysis allowed us to elucidate different contributions to 
the exchange couplings and understand the origin of these contributions on the microscopic 
level. Despite practical importance and broad interest to the HM ferromagnetism in Cr0 2 , 
the problem was far from being fully understood. There are several reasons for it. 

First, there is no single microscopic mechanism, which is primarily responsible for the 
ferromagnetism of Cr0 2 . Our analysis clearly shows that it is a joint effect of several 
contributions, of very different origins, and besides conventional double exchange in the 
t 2g band, there are other magnetic interactions, which are equally important in stabilizing 
the ferromagnetic ground state of Cr0 2 . They include direct exchange interactions, the 
interactions between t 2g and magnetically polarized oxygen 2 p band, as well as higher order 
effects in the (AE)^ 1 expansion for the magnetic energy of the t 2g band. 
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Second, the description of interatomic exchange interactions in Cr() 2 may have many 
traps, because the behavior of these interactions strongly depends on the method in use, 
which may lead to different conclusions. Particularly, if one sticks to static methods, which 
totally neglect the effect of correlation interactions on the magnetic properties (such as unre¬ 
stricted HF approximation), the solution of the problem may look very easy and the robust 
HM ferromagnetism emerges already in the minimal model, consisting of the f 2fl bands. 
However, this “easy solution” appears to be largely incomplete, as it becomes clear after 
considering the correlation interactions. Moreover, one should be most careful with the use 
of additional approximations for treating the correlation interactions, because some of these 
approximations may lead to unphysical results. For instance, by using the random-phase 
approximation for the correlation energy and treating this problem in the spirit of the OEP 
method with some static local potential, one can easily obtain an insulating solution, which 
suppresses the tendencies towards ferromagnetism. This curious example also demonstrates 
the importance of dynamic correlations. 


The most reliable technique for dealing with this kind of problem is the dynamical mean- 
field theory. In the present work, we have employed the new realization of this method, 
which is based on the exact diagonalization solution of the quantum impurity problem, 
performed ‘on-the-fly’. The use of this numerically advanced algorithm enabled us not 
only no to solve the standard DMFT equations, but also to study in many details the 
behavior of interatomic exchange interactions in the frameworks of this method. Our study 
provides an important insight into the origin of HM ferromagnetism in Cr0 2 . It clearly shows 
that, besides conventional processes, related to the change of the kinetic energy of electrons 
in the t 2g band, the realistic microscopic model for Cr0 2 should also include the direct 
exchange interactions and the magnetic polarization of the oxygen 2 p band. Finally, we 
have proposed how the latter two contributions can be evaluated using results of electronic 
structure calculations in the local-spin-density approximation. Thus, our work provides the 
firm microscopic basis for understanding the magnetic properties of Cr0 2 - the canonical 
and technologically important half-metallic ferromagnet. 
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